
	// Table A8, columns 3-4
	use chns_food, clear
	merge m:1 month using month_var_pca
	drop if _merge==2
	drop _merge
	keep if urban==0
	summ d3kcal
	egen std_food = std(food)
	reghdfe d3kcal Q4  , absorb(t1 wave urban hhid) cluster(hhid)
	reghdfe d3kcal std_food  , absorb(t1 wave urban hhid) cluster(hhid)


	use ceps_breastfed, clear
	
	// Table A21
	sum breast breast_4m junior water hhincgross std_hh_asset if  male==0 
	sum breast breast_4m junior water hhincgross std_hh_asset if  male==1 
	
	// Table 10
	summ breast   if  male==0 
 	summ breast   if  male==1 
	summ breast   if  male==0 & junior==0
 	summ breast   if  male==1 & junior==0
	summ breast   if  male==0 & junior==1
 	summ breast   if  male==1 & junior==1

	reghdfe breast Q2 Q3 Q4   if  male==0 , absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast Q2 Q3 Q4   if  male==0 & junior==0, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & junior==0 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast Q2 Q3 Q4   if  male==0 & junior==1, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & junior==1 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	
	
	// Table A22
	summ breast_4m   if  male==0 
 	summ breast_4m   if  male==1 
	summ breast_4m   if  male==0 & junior==0
 	summ breast_4m   if  male==1 & junior==0
	summ breast_4m   if  male==0 & junior==1
 	summ breast_4m   if  male==1 & junior==1

	reghdfe breast_4m Q2 Q3 Q4   if  male==0 , absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast_4m Q2 Q3 Q4   if  male==1 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast_4m Q2 Q3 Q4   if  male==0 & junior==0, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast_4m Q2 Q3 Q4   if  male==1 & junior==0 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast_4m Q2 Q3 Q4   if  male==0 & junior==1, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast_4m Q2 Q3 Q4   if  male==1 & junior==1 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	
	// Table A23 Panel A	
	summ breast   if  male==0 & water==0
 	summ breast   if  male==1 & water==0
	summ breast   if  male==0 & water==1
 	summ breast   if  male==1 & water==1

	reghdfe breast Q2 Q3 Q4   if  male==0 & water==0, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & water==0 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast Q2 Q3 Q4   if  male==0 & water==1, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & water==1 , absorb( ybirth wave prov ) cluster(prov_ybirth)
	
	// Table A23 Panel B
	summ breast   if  male==0 & hhincgross<=med_hhincgross
 	summ breast   if  male==1 & hhincgross<=med_hhincgross
	summ breast   if  male==0 & hhincgross>med_hhincgross
 	summ breast   if  male==1 & hhincgross>med_hhincgross

	reghdfe breast Q2 Q3 Q4   if  male==0 & hhincgross<=med_hhincgross, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & hhincgross<=med_hhincgross, absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast Q2 Q3 Q4   if  male==0 & hhincgross>med_hhincgross, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & hhincgross>med_hhincgross , absorb( ybirth wave prov ) cluster(prov_ybirth)

	// Table A23 Panel C
	summ breast   if  male==0 &  hh_asset<=med_hh_asset
 	summ breast   if  male==1 &  hh_asset<=med_hh_asset
	summ breast   if  male==0 &  hh_asset>med_hh_asset
 	summ breast   if  male==1 &  hh_asset>med_hh_asset

	reghdfe breast Q2 Q3 Q4   if  male==0 & hh_asset<=med_hh_asset, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & hh_asset<=med_hh_asset, absorb( ybirth wave prov ) cluster(prov_ybirth)
	reghdfe breast Q2 Q3 Q4   if  male==0 & hh_asset>med_hh_asset, absorb( ybirth wave prov ) cluster(prov_ybirth)
 	reghdfe breast Q2 Q3 Q4   if  male==1 & hh_asset>med_hh_asset , absorb( ybirth wave prov ) cluster(prov_ybirth)
	